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The origin of uniaxial and hydrostatic pressure effects on T c in the single-layered cuprate su- 
perconductors is theoretically explored. A two-orbital model, derived from first principles and 
analyzed with the fluctuation exchange approximation gives axial-dependent pressure coefficients, 
dT c /dP a > 0, dT c /dP c < 0, with a hydrostatic response dT c /dP > for both La214 and Hgl201 
cuprates, in qualitative agreement with experiments. Physically, this is shown to come from a uni- 
fied picture in which higher T c is achieved with an "orbital distillation", namely, the less the d x 2_ y 2 
main band is hybridized with the d z 2 and 4s orbitals higher the T c . Some implications for obtaining 
higher T c materials are discussed. 

PACS numbers: 74.20.-z, 74.62.Bf, 74.72.-h 



I. INTRODUCTION 

In the physics of high-T c cuprates, optimizing their T c 
remains a fundamental yet still open problem. Empir- 
ically, important parameters that control T c have been 
identified for the cuprates, i.e., chemical composition, 
structural parameters, the number of layers, etc, besides 
the doping concentration. For the structural parameters 
specifically, several key quantities have been suggested: 
the bond length between copper and in-plane oxygen(Z, 
defined in FigQ} 1,2 , and the Cu-apical oxygen distance 

Now, the pressure effect is exceptionally valuable as 
an in situ way to probe the structure-dependence of T c . 
Regarding this, two general observations have been made 
for the cuprates: (i) T c tends to be enhanced under hy- 
drostatic pressure, while (ii) uniaxial pressures produce 
anisotropic responses of T c . More precisely, (i) T c has 
been shown to monotonically increase for pressure < 30 
GP a 12 : 13 . As for (ii), an a-axis compression generally 
raises T c {dT c /dP a > 0), while c-axis compression has 
an opposite effect (dT c /dP c < 0)12r— . Moreover, the 
magnitude of the pressure coefficient becomes smaller for 
materials having higher T c 's, as summarized in Fig. 3 of 
ref>^. The purpose of the present study is to theoreti- 
cally reveal the origin of these general trends, focusing on 
the single-layered cuprates for clarity, and to shed light 
on a possibility of further optimizing T c . 

Conventionally, the theoretical model primarily used 
for the cuprates is a one-band Hubbard model based on 
the d x 2_ y 2 orbital (or sometimes Cu-3d x 2_ y 2+0-2p a or- 
bital). Recently, we have show n 10 : 11 that the d z 2 or- 
bital component strongly mixes into the states on the 
Fermi surface in the relatively low-T c cuprates such as 



La2Cu04 (La214)i 7 - _ — , where the hybridization works 
destructively against d-wave superconductivity. While 
there have been some theoretical studies in the literature 
focusing on the role of the d z 2 orbital 7 -!^ 2 ^— , Refsj 10 : 11 
conclude that larger the level offset AE between the 
d x 2_ y 2 and d z 2 Wannier orbitals, higher the T c , where 
AE is governed by the apical-oxygen height and the inter- 
layer distance. One might then presume that the effects 
of uniaxial pressures can simply be captured in terms 
of the pressure-dependence of AE affected by the crys- 
tal field. However, we reveal in the present work that the 
physics is not so simple. We find that, while the variation 
of T c under pressure is indeed affected by AE, especially 
in the relatively low-T c cuprates, the large AE values in 
higher-T c cuprates such as HgBa2Cu04 (Hgl201) make 
their relevance to the T c variation smaller. We shall show 
that we have to turn our attention rather to the Cu4s 
level, which is raised with pressure, resulting in a less 
rounded (i.e., better nested) Fermi surface. This, along 
with the increase in the band width, is shown to cause 
a higher T c under pressure. These results can be unified 
into a picture in which higher T c can be achieved by the 
"distillation" of the main (i.e., d x 2_ y 2~) band, namely, the 
smaller the hybridization of other orbital components the 
better. 



II. FORMULATION 
A. Construction of the two-orbital model 

Our theoretical procedure is as follows. We first de- 
termine the lattice structure under uniaxial and hydro- 
static pressures from a first-principles band calculation 



with the Wien2k package 24 . From the band structure, we 
construct the maximally-localized Wannier orbital o 25 i 26 
to obtain the hopping integrals for a two-orbital tight- 
binding model that takes into account both the d x 2_ y 2 
and the d z 2 Wannier orbitals explicitly^. 



B. Many body analysis 



with the interaction matrices 
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In this two-orbital model, we consider the onsite intra- 
and inter-orbital electron-electron repulsive interactions, 
which are given, in the standard notation, as 
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where z, j denote the sites while /x, ^ the two orbitals, the 
electron-electron interactions comprise the intraorbital 
repulsion U, interorbital repulsion U' , and the Hund's 
coupling J(= pair-hopping interaction J' ). Here we take 
U = 3.0 eV, V = 2.4 eV, and J=0.3 eV ? . These values 
conform to widely accepted, first-principles estimations 
for the cuprates that U is 7-10t (with t ~ 0.45 eV), while 
J, J' ~ 0.1C/. Here we also observe the orbital SU(2) 
requirement, U' = U — 2 J. 

To study the superconductivity in this multi-orbital 
Hubbard model, we apply the fluctuation exchange 
approximation(FLEX)2i~— . In FLEX, we start with 
Dyson equation to obtain the renormalized Green's func- 
tion, which is, in the multi-orbital case, a matrix in the 
orbital representation as G^i 2 , where l\ and li are orbital 
indices. The bubble and ladder diagrams constructed 
from the renormalized Green's function are then summed 
to obtain the spin and charge susceptibilities, 



Xs(q) 



Ml) 



X°(q) 



l-S X °(q) 



x°(q) 



(2) 



(3) 



where q = (q, iu> n ) with wave vector q and with Matsub- 
ara frequency iuj n = (2n + Vjirk^T, and the irreducible 
susceptibility is 



A./i,Z 2 ,^3,^ 



(q) =Y,Ghh(k + q)G hh (k), 
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With these susceptibilities, the fluctuation-mediated ef- 
fective interactions are obtained, which are used to calcu- 
late the self-energy. Then the renormalized Green's func- 
tions are determined self-consistently from Dyson equa- 
tion. The Green's functions and the susceptibilities are 
used to obtain the spin-singlet pairing interaction in the 
form 



V s (q) = ls X s(q)S ~ l -C X c{q)C + \{S + C), (7) 



and this is used in the linearized Eliashberg equation, 
T 



AA »' ( ' i ;= "irE E V«iW'(8) 

X G hh (k - q)A hh (k - q)G hh (q - k). (8) 



The superconducting transition temperature, T c , corre- 
sponds to the temperature at which the maximum eigen- 
value A of the Eliashberg equation reaches unity, so that 
A at a fixed temperature can be used as a measure for 
T c . T c of the Hg cuprate is experimentally about three 
times higher than in La cuprate^, so we calculate A by 
putting T = 0.01 eV for La and T = 0.03 eV for Hg 
for a clearer comparison. As we shall see, the eigenval- 
ues discussed in the present study are away from unity 
(i.e., the temperature is higher than T c ) due to the lim- 
itation in the number of Matsubara frequencies and the 
fc-point meshes. Therefore, for the La cuprate in par- 
ticular, we restrict ourselves to qualitative argument for 
the T c variation under pressure. For the Hg cuprate, on 
the other hand, we can go down to lower temperatures 
(T ~ 0.01) where the eigenvalue approaches unity, and 
we have checked that the conclusions drawn from the 
T = 0.03 calculation hold also for T ~ 0.01. Moreover, 
we estimate dT c /dP for Hg with the T ~ 0.01 results, 
as will be discussed in the final part of the paper. We 
fix the total band filling (number of electrons/ site) at 
n = 2.85, for which the filling of the main band amounts 
to 0.85 (15 % hole doping). We take a 32 x 32 x 4 k- 
point meshes for the three-dimensional lattice with 1024 
Matsubara frequencies. 
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TABLE I: Structural and electronic parameters obtained from 
the first-principles(th) and experiments(exp) in ie i 32 ^ 33 . 





La(exp) 


La(th) 


Hg(exp) 


Hg(th) 


a [A] 


3.78 


3.76 


3.88 


3.84 


CO [A] 


13.2 


13.1 


9.51 


9.58 


ho [A] 


2.42 


2.41 


2.78 


2.81 


^La,Ba [A] 


1.85 


1.81 


1.92 


1.88 


V [A 3 ] 


189 


184 


143 


141 


AE [eV] 


0.857 


0.861 


2.16 


2.305 


' x 2 —y 2 


0.363 


0.357 


0.419 


0.411 


W [eV] 


4.14 


4.23 


4.06 


4.19 



III. CALCULATION RESULTS: 
PRESSURE 



UNIAXIAL 



A. Crystal structure under pressure 

Let us begin with the case of uniaxial pressure. We 
first vary the lattice constants and calculate the total en- 
ergy E to t ■ This is fit by the standard Burch-Marnaghan 
equation 31 to determine the most stable structure with 
a unit cell volume V = Vo, the a- axis lattice constant 
a = ao, and the c-axis c — cq. For simplicity, we retain 
the tetragonal symmetry throughout, i.e., b — a (so that 
the a-compression is actually biaxial). We show in Ta- 
ble I the lattice parameters, ao, cq, ho, ^La,Ba(La or Ba 
height measured from CuC>2 plane) and Vo, obtained for 
the La and Hg cuprates. The results are in good agree- 
ment with experimental values for the optimally doped 
compound o 32 i 33 . We then relax the structure perpendic- 
ular to the compression direction, namely, we allow the 
lattice constant in that direction to relax to obtain the 
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FIG. 1: Bottom left: schematic variation of the d z 2 orbital 
level with respect to that for d x 2_ y 2 under uniaxial pressure 
(top left inset). Bottom right: the shift of Cu4s level under 
hydrostatic pressure and its effect on the Fermi surface. 
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FIG. 2: For uniaxial compressions the eigenvalue A of the 
Eliashberg equation is plotted against a/ao (a) or c/co (b). 
Triangles (circles) indicate the result for the La (Hg) cuprates. 
Arrows depict the contributions (see text) to the A variation 
from AE, W, and r x 2_ y 2, respectively, at a/ao, c/ Co = 0.975. 
Lines are guide for the eye, with the dashed horizontal ones 
indicating the original values. 



value that gives the lowest energy. 

Figure [5] plots the eigenvalue of the Eliashberg equa- 
tion A against the lattice compression a/ao a n d c/cq for 
each compound. The result shows that (i) in both com- 
pounds A increases as a/ao is reduced, while decreases 
as c/cq is reduced, and (ii) the absolute value of the 
variations of A is larger in La than in Hg. These fea- 
tures are in qualitative agreement with the experimental 
observations summarized in Fig. 3 of ref J^, which shows 
8T c /dP a > and dT c /dP c < for both materials, with 
larger \dT c /dPi\ in La than in H g 14 ' 15 . To be more pre- 
cise, while the compressibility in the a direction is nearly 
the same between the two materials, that in the c di- 
rection is about three times larger in Hg than in L a 34 i 35 
(dc/dP c | H g — 3dc/dP c | La), but even if we take this into 
account, we find that dX/dP c is still larger for La than 
for Hg in our calculation. 



B. Contribution from the d z 2 orbital: AE 

Now we want to pinpoint the origin of this T c varia- 
tion against uniaxial pressures. In both materials, AE = 
Edx 2 ~y 2 — Edz 2 increases as a/ao is reduced, while it de- 
creases when c/ Co is reduced. This is natural, since the a 
(c) reduction pushes the in-plane (out-of-plane) ligands 
toward Cu, resulting in a larger (smaller) crystal-field 



splitting and hence larger (smaller) AE^, as schemati- 
cally depicted in Fig[T] One might then expect that this 
alone is the origin of the T c variation, since AE and T c 
are positively correlated^ .. To see if this is indeed the 
case, we have considered a case where we increase AE 
alone to its value at a/aa = 0.975 or c/cq = 0.975, and 
obtain A with FLEX. The results are indicated in Fig[2] 
with arrows labeled as "AE" . In La, the resulting A is 
very close to those obtained for the actual compression, 
which implies that the main origin for A variation un- 
der uniaxial pressure comes from AE. By contrast, for 
Hg, the AE contribution is too small to account for the 
actual A variance (see the blowups in Figj2]). 



C. Contribution from the 4s orbital: r x i_ y i 

The reason for this is that in Hg, AE is ~ 2.5 times 
larger than in La (table H]) , so that the effect of the d z i 
orbital is tiny, while the contribution to the T c varia- 
tion coming from other changes in the electronic struc- 
ture become comparable with that from AE, In par- 
ticular, we focus on the change in the energy difference 
AE S between Cu4s and Cudj;2_ y 2 orbitals. In fact, it 
has been shown that the Cu4s orbital, which is implic- 
itly included in the d x 2_ y 2 Wannier orbital in the present 
scheme, affects the second (£2) and third (£3) neighbor 
hopping a 4 ' 6 ' 10 ! 11 . Note that the 4s orbital can be inte- 
grated out (implicitly included in the Wannier orbitals) 
prior to the many-body analysis, since the 4s orbital sits 
in energy well away from the Fermi level in contrast to 
the d z 2 orbital (FigJTJ- 10,11 . Smaller AE S results in larger 
r x 2_ y 2 = (|£ 2 | + |£3|)/|£i| within the d x 2_ y 2 orbital sec- 
tor, resulting in a more rounded Fermi surface, which 
degrades d-wave superconductivit y 10 ' 36- — . 

To show how the roundness varies with AE S , we con- 
sider a three-orbital model which explicitly includes the 
Cu 4s Wannier orbitals for the Hg cuprat e 10 ' 11 , and 
show in FigJ3] the Fermi surface for various values of 
AE S = Equ4s — Ecu3d !c 2_ 2- We stress here that, while 
larger AE and larger r x 2_ y 2 (or smaller AE S ) both give 
more rounded Fermi surface, their effects on T c are op- 
posite. Under pressure, AE S is enhanced, which in turn 
reduces r x 2_ y 2. In Fig[2J we show the effect of hypothet- 
ically reducing r x 2_ y 2 down to its values at a/ao = 0.975 
or c/co = 0.975. While the effect of r x 2_ y 2 is much 
smaller than that of AE in La, the two effects are found 
to be comparable in Hg. 



D. Contribution from the band width: W 

In addition to AE and r x 2 _ y 2 , the band width W (the 
energy difference between k — (0,0) and (71", it)) of the 
main band is also altered by pressure. In La the change 
in A due to the modification of W is small compared to 
that arising from AE, but in Hg the W contribution is 
comparable with those from AE and r x 2_ y 2, which in 



fact provides a full understanding of the net A variation 
under uniaxial pressure. Namely, the a (c) reduction re- 
sults in an increase (decrease) of the band width as ex- 
pected, which enhances (suppresses) T c . The increase of 
the band width results in a suppression of U/W, hence 
the electron correlation effect. This reduces the pairing 
interaction, while the self-energy correction due to the 
spin fluctuations is reduced at the same time. The former 
has an effect of enhancing T c , while the latter suppresses 
superconductivity. In the case of Hg compound, the lat- 
ter effect supersedes the former, resulting in an enhanced 
T c . 

It should be noted that the contribution from r x 2_ y 2, 
while relatively small for uniaxial compression, enhances 
T c for both of the a- and c- axis compressions in marked 
contrast with the contributions from AE and W. This 
will become important in our analysis for hydrostatic 
pressures below. 



IV. CALCULATION RESULTS: HYDROSTATIC 
PRESSURE 

A. La2Cu04 

Having identified the ingredients that determine the T c 
variation against uniaxial pressures, let us now move on 
to hydrostatic compression. Here we optimize the lattice 
structure at a fixed unit cell volume V(< Vq) by vary- 
ing Poisson's ratio, which we fit to the Burch-Marnaghan 
equation to obtain the most stable E tot . Notably enough, 
for hydrostatic pressures A in Fig|4] increases with the 
volume compression in both materials. This result qual- 
itatively agrees with experimental result s 12 ! 13 . To un- 
derstand its mechanism, we can, as done above for uni- 
axial pressures, decompose the pressure effect on A into 
the contributions from AE, W, and r x 2_ y 2 (arrows in 
Fig. Hfa)(b)). We can then realize that the variation 
of AE against hydrostatic pressure is not as straight- 
forward as in uniaxial pressures. Namely, we can look 
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FIG. 3: The Fermi surface of the three-orbital model of the 
Hg cuprate for values of AE S hypothetical varied from 6eV 
(nearly original value) to 12eV. 
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FIG. 4: For hydrostatic pressure applied to La(Hg) cuprates in the left(right) column, (a),(b): The eigenvalue A of the 
Eliashberg equation plotted against the volume compression V/Vq. Arrows are as in Fig[2]for V/Vo = 0.90. (c),(d): the value 
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Lines are guide for the eye. 
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at AE along with the "aspect ratio" ho/l against the 
volume reduction in Fig|4jc)(d), where ho is the apical 
oxygen height and I the in-plane Cu-0 distance. Under 
hydrostatic pressure, ho/l decreases in both materials 
because of the larger compressibility along the c direc- 
tion. One might then expect that this would reduce the 
crystal field splitting and hence AE, but actually this is 
by no means always the case. In fact, AE increases with 
pressure for La, which is because the Cu-0 distance de- 
creases, resulting in a larger crystal-field effect. Thus the 
T c enhancement in La mainly comes from the increase of 
AE. 



B. HgBa 2 Cu0 4 

The above argument for La does not directly apply to 
Hg, since the original apical-oxygen height is larger, so 
that there is more room for the CuO octahedron to shrink 
along the c-axis than in La. Therefore, the ho/l reduc- 
tion is larger, resulting in a nearly constant AE against 
the decrease of V/Vq. This further makes AE irrelevant 



to the T c variation in Hg under hydrostatic pressures. 
As seen in Fig|4] with arrows, main contributions to the 
T c enhancement come from W and r x 2_ y 2, with simi- 
lar magnitudes. As shown in FigJ3Je)(f), the decrease of 
r x i_ y i originally comes from an increase of the level off- 
set AE S introduced above. The relatively large enhance- 
ment of the Cu4s level under hydrostatic pressure can 
be understood from Fig IT] (right), where all the ligands 
approaching Cu push up the energy level of the extended 
and isotropic Cu4s orbital to a larger extent than for the 
localized and anisotropic Cu3<i orbitals. Thus a message 
here is the hydrostatic and uniaxial pressures exert sig- 
nificantly different effects. Specifically, the importance 
of r x 2_ y 2 becomes prominent in Hg in hydrostatic pres- 
sure because the r x 2_ y 2 -contribution is positive for both 
of a- and c-axis compressions, while VF-contribution has 
opposite effects as shown in Figj2j 

As for the band- width effect, we have found here that 
Hg exhibits an effect opposite to La for the present 
electron-electron interaction strength. To elaborate this, 
we have performed a FLEX calculation for various in- 
teraction strengths over 6 < U/t < 10, and found that 



increasing the band width always results in an enhanced 
A in Hg within the considered compression range, while 
in La a similar effect is obtained only for 8 < (U/t), 
with the effect reversing for smaller U. We have further 
noticed that this "sign change" in the band-width effect 
against U is peculiar to the systems having smaller AE. 
At any rate, the band-width effect is much smaller than 
the effect of AE in La, so that the effect of pressure- 
dependence of U does not affect the present conclusion. 



C. Order of magnitude of dT c /dP 

Let us finally comment on the relation between the A 
variation for hydrostatic pressures and the T c enhance- 
ment in the actual pressure experiments. To see this, we 
have extended our calculation to lower temperatures for 
Hg, where A becomes closer to unity (i.e., T approaches 
T c ). We find A ~ 0.86 at T = 0.01 eV for V = 0.9V , 
and the same value of A attained at T = 0.0088 eV for 
V = Vq, so the temperature difference (a rough estimate 
of AT C ) amounts to ~ 14 K. Since the compressibility is 
- 0.01 GPa" 1 ^, this implies dT c /dP ~ 1 K/GPa, which 
has the same order of magnitude found experimentally 14 . 



tions are revealed to come instead from the Fermi sur- 
face roundness governed by the Cu4s orbital as well as 
the variation of the band width W. These effects coming 
from the electronic structure in the multi-orbital systems 
can be unified into a single picture in which the orbital 
distillation of the main band results in a higher T c . 

The present study can also shed light on a materials- 
science avenue for optimizing T c . The strategy for en- 
hancing T c , as conceived here, is: (1) keep the level offset 
between the d x 2_ y 2 and d z 2 orbitals large (ideally, larger 
than U as shown in Fig [l] left), (2) expand the level offset 
between the Cu4s and the Cu3d x 2_ y 2 as much as possi- 
ble — this makes the Fermi surface more nested (figlU 
right), and (3) tune the band width W to a moderate 
value. In this sense, it is important to keep the distance 
ho between apical oxygen and Cu atom, and it is also 
important to decrease the in-plane Cu-0 bond length 
I. In other words, the desired situation for optimizing 
T c should have an a — b biaxial chemical pressure which 
reduces the length I from those in existing compounds, 
with the value of ho kept high. This may be coupled to 
the possibility of the level offset AE S controlled indepen- 
dently of AE by tuning length l la . 
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